Behavior of the thermopower in amorphous materials at the metal-insulator transition 
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We study the behavior of the thermal transport properties in three-dimensional disordered sys- 
tems close to the metal-insulator transition within linear response. Using a suitable form for the 
energy-dependent conductivity a, we show that the value of the dynamical scaling exponent for 
noninteracting disordered systems such as the Anderson model of localization can be reproduced. 
Furthermore, the values of the thermopower S have the right order of magnitude close to the transi- 
tion as compared to the experimental results. A sign change in the thermoelectric power S — as is 
often observed in experiments — can also be modeled within the linear response formulation using 
modified experimental a data as input. 
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I. INTRODUCTION 

Transport phenomena in disordered quantum systems 
have been studied for many years,EH3 yet many open prob- 
lems remain. One focus of these investigations is the 
metal-insulator transition (MIT). This quantum phase 
transition from a good conducting material to an insu- 
lator may be induced by disorder due to localization^ or 
by interactions such as electron-electron interactions and 
electron-lattice couplingfltj In three-dimensional (3D) 
amorphous materials the MIT is mainly attributed to 
disorder .EJ For example, in heavily doped semiconductors 
the disorder is brought about by the random distribution 
of dopant atoms in the crystalline host. However, indi- 
cations of electron-electron interactions have also been 
found, e.g., in the d.c. conductivity a (or resistivity 
p = 1/cr) in doped semiconductors in both metallidj and 
insulating regimes. □ 

A further open problem is the behavior of the ther- 
moelectric power S or the Seebeck coefficient of disor- 
dered materials near the MIT. In amorphous alloys and 
both compensated Si:(P,B) and uncompensated Si:P, S 
continuously changes from negative to positive values or 
vice versa at low temperature T. This corresponds to a 
change of thermal conductors from electrons to holes or 
conversely and has been attributed to electron-phonon 
interaction in amorphous alloysfffl On the other hand 
in heavily doped semiconductors the sign change is be- 
lieved to be caused by electron-electron interactions or 
attributed to the existence of local raaaaetic moments 
and their interactions with electrons BBB This conclu- 
sion is based on the suppression of the anomalous be- 
havior by a magnetic field. ua We remark that the sign 
change in S is also qbaecyed in metals, high-T c materi- 
als and quasicrystals£2rO Analytical treatments of met- 
als as a degenerate free-electron gaSptaking into account 
inelastic scattering with phononsEj lLj as well as numer- 
ical considerations-incorporating electronic correlations 
in superconductorsE3 have also been shown to generate a 



sign change in S. But in these systems the sign change 
occurs at a T value which is 2 orders of magnitude higher 
than in disordered systems. Note that S is mainly due 
to two distinct effects: (i) the diffusion of the charge car- 
riers and (ii), the net momentum transfer from phonons 
to carriers. Ell But for T < 0.3 K as considered in this 
work, the diffusive part of the thermoaower dominates 
that of the phonon-drag contributionO'Eij Hence, from 
this point on in this paper S denotes only the diffusion 
thermopower. 

The prototype for a theoretical description of 3D dis=. 
ordered systems is the Anderson model of localization.EB 
Near the MIT at T = 0, a behaves as 




\E V \ < E c 
\E F \ > E c 



(1) 



where E-p is the Fermi energy, E c is the mobility edge 
which separates the extended conducting states from lo- 
calized insulating states, and v is a universal critical 
exponent a By using Eq. ([!]) for a in a linear response 
formulation the behavioij— of— the thermoelectric trans- 
port properties such as S'jHj'Ea thc-tjicrmal conductivity 
lffl~Eland the Lorenz number L$3Q at the MIT have 
been computed. Moreover, similar to cr, the quantities 
S, K and Lq have also been found to obey scaling.O The 
scaling form of the dyaaquffi] conductivity a close to the 
MIT in 3D is given asm 



<r(t,T) 

rpl/x 



t 



(2) 



Here t measures a dimensionless distance from the crit- 
ical point, such as t — (Ep — E c )/E c , the correlation- 
length exponent v in 3D is equivalent to the conductiv- 
ity exponent as given in Eq. (!]), and z is the dynamical 
exponent fl For a noninteracting system such as the An- 
derson model, one expects z = d in d dimcnsions.a But, 
instead of obtaining z = 3 in the scaling form of a, one 
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finds zv 



i E3SE3 



In addition to this discrepancy 



turns out to be at least one order of magnitude largerc 
than the experimental results in doped semiconductors^ 
and in amorphous alloys .EIIj Furthermore, the sign change 
in S cannot be explained using the Anderson model and 
Eq. (Q). One may argue that the discrepancies between 
the transport calculations and the experimental measure- 
ments are due to the absence of interactions in the An- 
derson model. Indeed, interactions may influence the be- 
havior of the thermoelectric transport properties. Yet we 
emphasize that the neglect of interactions in the Ander- 
son model is not entirely inconsistent with the experimen- 
tal situation in 3D amorphous materials. For example, 
recent measurements in Si:P yield a scaling with z ss 3 
and v w lr3 ThiS|-agrees with z = d as predicted by the 
scaling argumentsonJ for noninteracting systems. 

The goal of this paper is to show that the correct value 
of z, the right order of magnitude of S at the MIT, and 
perhaps even the sign change in S at low T, can be de- 
scribed within a linear response formulation using the 
noninteracting Anderson model of localization. However, 
in order to do so, we have to use a more suitably chosen 
energy-dependent er c instead of Eq. (Q) . After a brief re- 
view of linear transport theory, we construct a new form 
for er c as a function of energy E and T from experimental 
data. By using this model data as input for the linear 
response formulation, we compute the temperature de- 
pendence of S, K, Lo and also a and show that they 
have the expected qualitative and quantitative behavior 
close to the MIT. Finally, we show that a small variation 
in <r c (E,T) can change the sign of S. This effect cannot 
be produced simply by varying the density of states g or 
the chemical potential fi(T). 



II. LINEAR THERMOELECTRIC TRANSPORT 
THEORY 



Equation (H) then yields the thermoelectric power S 
which relates VT to E, 



S 



\e\TLu 



(6) 



The sign of S determines whether the thermal carriers 
are electrons or holes. Using the Som fiM^ Jd] expansion 
for \E F - E c \> k B T, S is given asH 2 ^ 



7r 2 fc|T dlncr(£0 



3|e| 



dE 



(7) 



E — Ep 



where E-p is the Fermi energy, &b is Boltzmann's constant 
and a(E) is assumed to be a slowly varying function on 
the scale-t>f &bT. Equation (Q) is also known as the Mott 
formulaS 

The thermal conductivity K determines the contribu- 
tion to (J2) stemming from VT. Using Eqs. (||) and 
(J2) we obtain K in terms of the kinetic coefficient ad 



K = 



L22L11 — L<i\L\<2. 
\e\ 2 TL n ' 



(8) 



From- 
thatE 



le definition of the Lorenz number it follows 



Ln 



L22L11 — L21L12 



(9) 



In metals at room T, Lq 



7r 2 /3.0. It also takes on the 
2re the electrons 



same value at T < 10 K in metals wl 
suffer no inelastic scattering processes.EH 

The primary consideration then in determining a, S, 
K and L is to calculate Ly. Under the assumptions 
that the system is noninteracting and inelastic scatter- 
ing processes are absent, L t j are given-in— the Chester- 



Thellung-Kubo-Greenwood formulationo~E3 as 



In the presence of a small temperature gradient VT, 
the electric current density (j 1 ) and the heat current den- 
sity (J2} induced in a system are given (to linear order) 
as 



Hi) = |e|-* (\e\LaE - L a T- l VT) 



(3) 



where e is the electron charge and E is the induced elec- 
tric field. Lij are the kinetic coefficients. Since we do not 
consider the presence of a magnetp^field in this work, the 
Onsager relation 



Lji holdsrj Ohm's law, 



implies that in Eq. 



Oi) = <rE, 



a = Lu 



(4) 



(5) 



The flow of thermal conductors due to VT is counter- 
acted by an electric force arising from E making (ji) = 0. 



Ui= \ A{E)[E - n(T)] 



i+j-2 



dE 



dE, 
(10) 



where i,j = 1,2, /i(T) is the chemical potential, 
f(E,fi,T) is the Fermi distribution function and A(E) 
contains all the system-dependent features. 

Lastly, we note that the T dependence of /j, can be 
obtained for noninteracting systems from 



n((i,T) 



dEg(E)f(E, f i.T) 



11) 



where n is the number density of electrons and g is the 
density of states.Eil Knowing g and keeping n constant, 
we find numerically that n(T) ~ T 2 in the 3D Ander- 
son modfiL with an increased effective mass due to the 
disorder .E3 
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III. A PHENOMENOLOGICAL APPROACH 

There are only two parameters that are model depen- 
dent in the transport theory discussed in Sec. |j| These 
are A(E) and fi(T). In order to determine the behav- 
ior of the thermoelectric transport properties close to 
tho-^aderson MIT, A(E) in Eq. (jig) has usually been 



equal to the critical behavior of a given by 
Eq. (0). As mentioned in the introduction, this leads to 
the unphysical value for z — l/v and therefore an un- 
physical frequency and T dependence of a. The main 
reason for this behavior is easily understood: there is no 
T dependence in Eq. ([!]) and thus all T dependence in 
Eq. ( 3) is due to the broadening of the Fermi function in 
Eq. (1C|) with increasing T. Thus in order to model the 
correct T dependence, we should add to Eq. ([!]), valid 
at T = 0, the desired T dependencies such as cr oc T x l z 
in the metallic and a oc exp (— T) in the insulating (say, 
variable-range-hopping) regimes. Such a purely theoreti- 
cal model for a c (E, T) will then incorporate a multitude 
of constants that can be adjusted to fit the experimental 
results. Of course this is of limited practical use since the 
validity of these fitting parameters is hard to justify. 

Here we will instead use as input for cr c {E 1 T) recent ex- 
perimental data obtained by Waffcnschmidt et aZ.ES who 
measured a in Si:P at the MIT under uniaxial stress. 
Their data yield good scaling of a according to Eq. (||) 
with a dynamical exponent z = 2.94±0.3 and v =p..LjL0.1. 
These values agree with the scaling argupieptBoEa and 
reasonably well with the numerical resultsEa for non- 
interacting systems. 

In the a(t, T) scaling of Ref. t = (s — s c )/s c where 
s is the stress and s c the corresponding value at the tran- 
sition. We sample their scaled data for several values of 
(t, T) and fit a spline curveEj a c to these points in order 
to get a smooth functional form for a(t,T). Transform- 
ing the spline a c as a function not only of T but also of E, 
we set t = (E - E c )/E c . Finally, we substitute a c (E,T) 
for A(E) in Eq. (|10|) . Consequently, the thermoelectric 
transport properties are directly obtained from Eqs. (|^)- 

In this paper we consider temperatures from 0.01 K 
to 0.2 K. Far from the transition we could not probe 
lower than T < 0.02 K. This is due to the limited input 
data and consequently a limited range of the spline func- 
tion that generated a c (E,T). The unit of a c is taken 
as f2 _1 cm _1 consistent with the experiments. The E 
scale is (arbitrarily) fixed at 1 meV which is the order 
of magnitude of the binding energv of an isolated donor 
in a heavily doped semiconductor .3 In order to .compare 
with the previous results in the Anderson modelEj we let 
E c = 7.5. We emphasize that this value is of no sig- 
nificance and can be assigned (nearly) arbitrarily. The 
important point to consider is the location of the Fermi 
energy Ep with respect to E c . This distinguishes the 
electronic regimes. Thus, the metallic, critical and insu- 
lating regimes are identified as \Ep\ < E c , Ep = E c , and 



\Ep\ > E c , respectively. Usually, fi(T) is derived in Eq. 
( |lT| ) from g of the 3D Anderson model of localization. In 
the next section we shall also show the effect of using a 
different functional form of fJ-(T). 



IV. RESULTS AND DISCUSSIONS 

A. Temperature dependence of the thermoelectric 
transport properties 

Consistent with the dynamics of the experiment in Ref. 
p8| , we expect <r(T) ~ T 1 ! 7 " at the critical regime with 
z w 3. This is indeed the behavior of <r(T) close to E c 
as we show in Fig. [l| For \Ep — E c \ < 0.2 meV we 
obtain z — 3.2 ± 0.3. Note that <r(T) = Ly\ has been 
integrated according to Eq. ( |l0| ) over the energy range 
where df/dE > 10~ 20 meV. Thus our numerical calcu- 
lation of a is consistent since it reproduces closely the 
original result in Ref. |28|. If we plot the results in Fig. 
|l| with respect to (jj, — E c ) j 'E C T Y I VZ we obtain a rough 
scaling of a similar to Fig. 3 of Ref. |2§| . 

We next turn our attention to the thermoelectric power 
S. In the 3D Anderson model of localization, we know 
that wl 
regimcE ^ 2 1 



ing Eq. (|1|) one obtains S — > in the metallic 
ic3 while in the insulating regime -Si-dpes not 
approach zero but seems to diverge as T — > O.BI 2 J r ALJ 
MIT S is a constants of the order of 100 ^V/K.BlE 
In Fig. ^j we show that in the present approach S in the 
vicinity of the MIT is two orders of magnitude smaller 
compared to previous results for the Anderson model. 
The magnitude of S is in fact comparable to the exper- 
imental results in disordered systems.tffitl Furthermore, 
S — » as T — > in the metallic, critical and insulating 
cases. This behavior of S(T) in all electronic regimes 



was observedB'cil in amorphous Au x Sbi_ x and in amor- 
phous Gei_2;Au x . As indicated by the different lines in 
Fig. 2, S is in good agreement with Eq. (0) since a in 
Fig. 1 is smooth across the transition at finite T. Note 
that in order to evaluate Eq. (0) properly for the system 
considered here, the E dependence of the input spline 
a c (E,T) was used instead of a(T) from Fig. |[ We em- 
phasize that it is no contradiction that S is positive here 
but mainly negative in the doped semiconductors in all 
electronic regimes. In the energy regions close to E c > 
the charge carriers are holes instead of electrons as shown 
in Ref. |23|. S would be negative if we_had chosen the left 
mobility edge E c < for low fillingOo 

The corresponding T dependence of K is shown in Fig. 
||. We find that K — > as T — + in all electronic .cegimes. 
This is also the behavior of K using a c in Eq. (0) .El In the 
metallic regime K is larger than in the insulating regime 
since there are more heat carriers in the former. From the 
results of a and K in Figs. [I] and|| respectively, we obtain 
io- As shown in Fig. [|, Lq — > 7r 2 /3 as T — > whether 
it be in the metallic, critical or insulating regime. This 
is different from the results using Eq. (jy) for a. There 
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one obtains an Lq which depends on the conductivity 
exponent in the critical and insulating regimes while it 
approaches the universal value&il 7r 2 /3 only in the metallic 
regime SB Here we see no markedly distinct behavior in 
the metallic regime compared to the insulating regime. 
For \Ep — E c \ = 0.1, 0.2 meV, Lq in the metallic regime is 
less than its corresponding value in the insulating regime. 
For \Ep — E c \ — 0.5, 1.0 meV, Lq in the metallic regime 
is larger than its corresponding value in the insulating 
regime. 

In the calculation of cr(T), S{T), K{T) and L Q (T), we 
used a phenomenological construction of a c (E,T). Fur- 
thermore, we have assumed that the density of states g 
is the same as that of the 3D Anderson model of local- 
ization given in Ref. |2^. Since this g{E) is a smooth and 
(restricted to E > 0) monotonic function, /x(T) obtained 
from Eq. ( |ll|) is also smoothly and monotonically varying 
with T as described in Sec. O. 



B. Effects of a structured g 

We now consider the effects of a possible structure in 
g. We shall assume here that this structure corresponds 
only to variations in /i(T) and not inc. In Fig. |^ we show 
two examples of a modified jtx(T) in the critical regime. 
Example A has a pronounced maximum, while example 
B has both a maximum and a minimum. The height of 
the maximum in both examples A and B is ~ 0.1 meV. 
Note that this is significantly larger than the halfwidth of 
the bump which is < 0.005 meV. This is also true for the 
depth of the minimum in example B. Thus a small change 
in T corresponds to a large change in u(T). Applying 
these forms of /i(T) together with Eq. (jl]) for a repro- 
duces the same structures in S. For example, using form 
B of n{T) we obtain an S having both a maximum and a 
minimum in the same T interval as n(T). But S is still of 
the order of 100/iV/K while the variations are only of the 
order of 10/xV/K and not large enough to cause a sign 
change in S. On the other hand using the phenomenolog- 
ical construction of a c (E, T) yields even smaller changes. 
We observe variations in L\\ and L\i of less than 10% 
from their unmodified values. Consequently, wc find neg- 
ligible changes in S. Figure [2] would appear unmodified. 
Hence we conclude that even a large change in the density 
of states g and thus also in ^(T) is not sufficient to cause 
the change of sign for S as observed in the experiments. 
However, this weak dependence of S on g and /z(T) at 
least justifies our use of the simple Anderson density of 
states in the present paper. 



C. Effects of a structured a(E) 

Let us now assume that for small T there are non- 
monotonicities in a c (E,T) — althaugh these have not 
been observed in the experiments. BQ'u Thus we consider 



the case when there is a sizeable change in a c (E, T) in 
the region close to E c for small T. The correspond- 
ing 'bumps' in <r c (E,T) are shown in Fig. || with dif- 
ferent peak heights and with half-widths < 1 meV. 
For simplicity they are essentially quadratic functions 
of E and have been generated such that they decay 
quickly as exp(— T 4 ) with increasing T . The height of 
the bumps is < 1 f7 _1 cm _1 which is at least an order 
of magnitude smalkt-than the values observed for a in 
the measurements.BEj The lowest temperature studied 
is T = 6 mK and we shall only consider metallic and 
insulating regions with \Ep — E c \ < 0.1 meV. 

Our results in Fig. |?| using u c {E, T) with and without 
bumps indicate that there are only small variations in 
the slope of log(cr) and a ~ T 1 / 3 remains valid within 
the accuracy of these estimations. We remark that the 
lowest measured temperature in Ref. ^8] is 15 mK. From 
Fig. ^ we see that the variations for T > 15 mK are much 
smaller than for T < 15 mK. Hence, these variations in 
a could not have been observed in the experiments. 

In Fig. U we show how the bumps affect S. Even with 
the very small bump 3, S changes sign in the critical 
regime as T — > 0. As the bump increases this change be- 
comes more pronounced. The temperature T$ = q at which 
the sign change occurs is 0.1 K for bump 3, Ts=o = 0.2 
K for bump 2, and Ts=q ~ 0.4 K for bump 1. These 
results for T$ = o are still about one order of magnitude 
less than in semiconductors!™ and two .Qjeders of magni- 
tude smaller than in amorphous alloys.ErLJ Of course, as 
shown in Fig. ^[ Ts=o shifts to higher values as the bump 
height increases. Nevertheless, the minimum value of S 
for T < Ts=o has the same order of magnitude as the cor- 
responding maximum value of S in Si:(P,B) and in Si:P. 
We emphasize that the value of Ts = o of course depends 
on the energy unit chosen and thus will vary for systems 
with different bandwidths, e.g., a larger band width will 
give rise to a larger value of Ts=o- We observe a simi- 
lar sign change in the metallic regime but the depth of 
the minimum is smaller than in the critical regime. The 
Mott formula (0) with <r(T) given in Fig. ^ can readily 
model this behavior since <r(T) remains slowly varying 
even if a c (E,T) has a bump near E c . In the insulating 
regime, S has a shallow maximum and drops back to zero 
as T — > 0. This is different from experiment ,u where S 
changes sign and neither has a maximum nor minimum 
in the insulating regime. 

Unlike in S, there is no dramatic change in K as can 
be seen in Fig. ||. We find only negligible variations at 
T < 15 mK. This should be expected since there has also 
been hardly any change in the slope of cr(T) except for 
T < 15 mK. However, the small increase in a at T < 15 
mK in the metallic regime together with the minimally 
modified K leads to a drastic change in Lq even in the 
case for the smallest bump. The increase and decrease in 
a leads to a maximum and minimum in Lq, respectively. 
However, Lq still approaches the universal value 7r 2 /3 for 
T — > as demonstrated in Fig. fi~0|. 
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V. CONCLUSIONS 



In this paper, we have shown that the anticipated 
value of the dynamical scaling exponent z 3 as well 
as the right order of magnitude for the thermopower 
S w 1/iV/K at the MIT can be obtained when taking into 
account the expected T dependencies in addition to the 
simple scaling behavior of Eq. (Q) . Our approach is phe- 
nomenological in the sense that we have refrained from 
using fitting parameters and have rather taken experi- 
mental data as input. Using this data, we can explain the 
large deviations from experimental results as reported in 
the theoretical studies of Refs. |20|- |24| , |3"o| . We have shown 
that our results for S agree with those predicted by the 
Mott formula (Q) since we have used a er slowly vary- 
ing on the scale of k^T near the MIT. We emphasize, 
however, that for a disordered system where interactions 
are negligible, we should still expect the Anderson-type 
transition as given in Eq_JIL.at T = 0. Consequently, 
S w 100/iV/K at the MITdTB should again be expected 
and one should observe a large increase of 151 at very low 
T. However, such temperatures appear presently inac- 
cessible by experiment. 

As a further challenge, we considered the sign changes 
observed in S at low T. We found that even large varia- 
tions in the chemical potential fi(T) do not lead to a sign 
change in S. On the other hand, a variation in the input 
a c (E,T) data can give rise to such a sign change in S, 
while at the same time resulting in only small changes in 
the conductivity a. Hence we have effectively modeled 
the underlying physical reasons for the sign change - 
which have been attributed to electron-electron interac- 
tions or to the existence of local-magnetic moments and 
their interactions with .djgctronslfflj or to inelastic scat- 
tering with phononsEjtJ — by simply changing the in- 
put a c (E,T). Regarding a possible test for the existence 
of such a structured a c (E, T) , we have shown that the 
T variation of Lq is much more sensitive to the bumps 
than a. Thus we have been able to describe the main 
features of the critical behavior of S(T) although it re- 
mains unclear what might cause bumps in u c {E, T) close 
to E c . A microscopic and possible system-dependent ap- 
proach to the problem may eventually account for these 
abrupt changes in a c . Of course, if many-particle inter- 
actions and electron-phonon coupling are important, we 
no longer expect the fcasibilits2-p£-thc Chester-Thcllung- 
Kubo-Greenwood formulatior£jL3 used here. 
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FIG. 1. Numerical calculations for the electrical conduc- 
tivity, according to Eq. (|^), as a function of temperature T. 
The filled symbols represent the metallic regime \Ef\ < E c , 
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represent the insulating regime \Ef\ > E c . 
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